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ABSTRACT 

The asymptotic capacity at low input powers of an average- 
power limited or an average- and peak-power limited discrete- 
time Poisson channel is considered. For a Poisson chan- 
nel whose dark current is zero or decays to zero linearly 
with its average input power £, capacity scales like £ log i 
for small £. For a Poisson channel whose dark current is 
a nonzero constant, capacity scales, to within a constant, 
like £ log log j for small £. 

Index Terms — Asymptotic, Capacity, Low SNR, Pois- 
son channel. 

1. INTRODUCTION 

We consider the discrete-time memoryless Poisson channel 
whose input x is in the set Rq of nonnegative reals and whose 
output y is in the set Zq of nonnegative integers. Conditional 
on the input X — x, the output Y has a Poisson distribution 
of mean (A + x) where A > is called the dark current. We 
denote the Poisson distribution of mean C by 7^{(-) so 

With this notation the channel law Vt^( |-) is given by 

W{y\x)^Vx+^{y), xeM+yeZ+. (1) 

This channel is often used to model pulse-amplitude mod- 
ulated optical communication with a direct-detection receiver 
n]. Here the input x is proportional to the product of the 
transmitted light intensity by the pulse duration; the output y 
models the number of photons arriving at the receiver during 
the pulse duration; and A models the average number of extra- 
neous counts that appear in y in addition to those associated 
with the illumination x. 

The average-power constrain^ on the input is 

E[X]<£, (2) 

'The word "power" here has the meaning "average number of photons 
transmitted per channel use." If we denote by P the standard "power" in 
physics, namely, energy per unit time (in watts), then the notation of "power" 
in this paper is really r\PT j tko, where t) is the detector's quantum efficiency, 
T is the pulse duration (in sec), and tku is the photon energy (in joules) at the 
operating frequency w (in rad/sec). 



where £ > is the maximum allowed average power 

The peak-power constraint on the input is that with prob- 
ability one 

X<A. (3) 

When no peak-power constraint is imposed, we write ^ = oo. 

No analytic expression for the capacity of the Poisson 
channel is known. In 1 1 1 Shamai showed that capacity- 
achieving input distributions are discrete whose numbers of 
mass points depend on £ and A. In ||2] [31 Lapidoth and 
Moser derived the asymptotic capacity of the Poisson chan- 
nel in the regime where both the average and peak powers 
tend to infinity with their ratio fixed. 

In the present paper, we seek the asymptotic capacity of 
the Poisson channel when the average input power tends to 
zero. The peak-power constraint, when considered, is held 
constant and hence does not tend to zero with the average 
power. We consider two different cases for the dark cur- 
rent A. The first case is when the dark current tends to zero 
proportionally with the average power. This corresponds to 
the wide-band regime where the pulse duration tends to zero. 
The second case is when the dark current is constant. This 
corresponds to the regime where the transmitter is weak. 

Our lower bounds on channel capacity are all based on 
binary inputs. In some cases we show that this is asymptoti- 
cally optimal. Our upper bounds are derived using the duality 
expression (see f4l and references therein). An efficient way 
to compute asymptotic capacities at low average input pow- 
ers is to compute the capacity per unit cost |5 |. However, we 
shall see that, apart from one case (Equation (|7]l), the capacity 
per unit cost does not exist, namely, the capacity tends to zero 
more slowly than linearly with the average power. 

Among the results in this paper, the special case of zero 
dark current has been derived independently in ||6l [T) • 

The rest of the paper is arranged as follows: in Section |2] 
we state the results of this paper; in Section [3] we prove the 
lower bounds; and in Section |4] we sketch the proofs for the 
upper bounds. 



2. RESULTS 

Let C{X,£,A) denote the capacity of the Poisson channel 
with dark current A under Constraints (|2]i and (O 

C{X,£,A) = snp I{X;Y) 

where the supremum is over all input distributions satisfy- 
ing © and ©. 

When A is proportional to £, the asymptotic capacity of 
the Poisson channel as £ | is given in the following propo- 
sition. Note that this also includes the case where the dark 
current is the constant zero. 

Proposition 1 (Dark Current Proportional to £). For any c > 

and A e (0, oo], 

C{(£,£,A) ^ 

lim ; = 1. 

Sio £ log i 

Recall that, for any a, (3 > 0, the sum of two indepen- 
dent random variables with the Poisson distributions Vai-) 
and 7^/3 (•) has the Poisson distribution Va+fsi')- Thus, we 
can produce any Poisson channel with nonzero dark current 
by adding noise to a Poisson channel with zero dark current. 
Consequently, 

C(0, A) > C{c£,£,A), c,£,A>0. 

Thus, to prove Proposition [T] we only need to show the fol- 
lowing two bounds: 

^^ C{c£ £ A) c>0,^e (0,oo], (4) 

£IQ £ log J 



— CiO,£,A) ^ Ae{0,oo]. 

£io flogi - ' ^ ' ^ 



(5) 



We shall prove (|4]i in Section 13.11 and shall sketch a proof 
for (0) in Sectionl4Tn 



Remark 1. The bound (|5]l can also be derived by noting that 
the capacity of the Poisson channel with zero dark current 
under an average-power constraint only is upper-bounded by 
the capacity of the pure-loss bosonic channel, and by using 
the explicit formula 



of the latter 



,{£)^{l + £)\og{l + £)-£\og£ (6) 



Remark 2. Because the pure-loss bosonic channel with co- 
herent input states and direct detection reduces to a Poisson 
channel, the lower bound (|4|l and the achievability of its left- 
hand side using binary signaling combine with ^ to show 
that the asymptotic (quantum-receiver) capacity of the pure- 
loss bosonic channel is achievable with binary modulation 
(on-off keying) and direct detection. 



For a Poisson channel with constant nonzero dark current, 
we have the following result. 

Proposition 2 (Constant Nonzero Dark Current). For any 

A > 0, 



and 



(7) 



'-<X^9^M^<^^^^<2. (8) 
2 flog log T £io flog log i 



The proof of (|7]) is a simple application of the formula for 
capacity per unit cost |5 , Theorem 2]. The proof of the lower 
bound in ^ is in Section |372l and a sketch of the proof of the 
upper bound in (|8]l is in Section|42] 



3. THE LOWER BOUNDS 

The achievability results in this section are obtained by choos- 
ing binary input distributions and then computing the mutual 
informations. We denote by Q'° the binary distribution 



X = 



0, w.p. (1 -_p), 
C, w.p. _p, 



(9) 



where ^ > 0, p G (0, 1). If we choose the parameters Q and p 
in such a way that Constraints (|2]) and (O are satisfied, then 

C{\£,A)>I{Q'°,W). (10) 



3.1. Dark Current Proportional to £ 

In this subsection we shall derive Inequality (|4]l. To this end, 
we compute the mutual information /(Q*", W) for input dis- 
tribution given by (|9]i: 

/(Qb, W) = H{Y) - H{Y\X) 

OO 

= - ^((1 - p)Vx{y) +pVx+dy)) 
■\og{{l ~p)Vx{y)+pVx+dy)) 



{l-p)Y,Vx{y) \ogVx{y) 



t/=0 



/o(A,C,p) + /i(A,C,p), (11) 



where in the last equaHty we defined 

/o(A,C,p)=-((l-p)e-'+pe-(^+^)) 

■log((l-p)e-^+pe-(^+«)) 

oo 

h{XX,p) " ^Y.i'^l -p)Vx{y) +pVx+dy)) 
■\og{{l -p)Vx{y)+pVx+dy)) 

CO 

+ {l-p)Y,'P>^iy)^ogVx{y) 

oc 

Note that in the above decomposition we took out the terms 
corresponding to y — in all three summations to form 
/o(A, C,p) and collected the remaining terms in Ii (A, C, p)- 
We lower-bound Io{XX,p) as 



/o(A, C,P) > - (1 - p)Xe-^ - p{X + C)e-(^+';) 
> -X-p{X + 0- 

We lower-bound /i (A, C, p) as 
/i(A,C,p) 

OO 

= - ^((1 - p)Vx{y) + pVx+dy)) 



(12) 



(l-p)^PA(y)log 



y=l 



Vxjy) 
Vx+dv) 



1-p Vx{y) 



\ogp + log 1 



p 'P\+dy) 



)) 



{1 -p)Y^V^{y) log — 



.-A A" 



fA+t;) (A+C)" 



y=i 

1~P Vx{y) 



logp + log 1 



p 'P\+dy) 



^ l-p T>.{v) 



(1 - p)c ^A(y) +(1 - p) log ^ E '^>^^y)y 

v=i ^ v=i 



>~Y.^{l-p)Vx{y)+pVx+dy))^oZP 



= (l-p)(l-e-^)+p(l-e-(^ + C)) 
oo 

^Y.{{i-p)Vx{y)+prx+dy)) 

y=i 



1-P Vxjy) 
p 'Px+dy) 

C 

X 



+ (l-p)(l-e-'^)C-(l-p)Alog(^l 
= ((1 - p){l - e-^)+p{l - e-^>^+<y)) log - 



(I-P)^ f- (^a(2/))^ 



-(l-p)^pA(y) 



^e^-p ^2 (y) 

A + C 

+ (1 - - e-^)C - (1 - p)Alog (^1 + ^ 
((1 -p){l- e-^) + p(l - e_(^3) log i 



—1 — 



<e-C 



P 



2 .2 



■e^(^l-e - (1 -p) (1 - e"-^) 

^ ? ! ' <i <A 



- - (1 - p)Aiog (^1 + ^ 

<1 <A 

> (1 - - e-^) log -+p{l- e-^) log i 
p p 

~~e^-X-XC~{\-p)X\og[l + i^. (13) 



A 



Choose any C, G (0, yl] and, for small enough £, let p = 
£/C. Then the distribution ^ satisfies both Constraints (|2j 
and Q. Let A = <£. Using (fT2l i we can bound the asymptotic 
behavior of /q ( A, C, p) as 



/o(cf,C,f) 
lim -r^ ^ — — > — lim 



.. f (eg + C) 



— lim ■ 



£10 flog I fioflogi fio £logi 



= 0. 
(14) 

Similarly, using (fTSl) we can bound the asymptotic behavior 
of /i(A, C,p) as 



lim 



h (cf,C,f) 



> 



£ia log J C 
Combining ( fTOl l. ( fTTT l. ( fT4l l. and ( fTSl l we obtain 



(15) 



lim ■ 

£10 £ log ; 



> 



c 



, forallC e (0,^]. (16) 



= 1-6- 



=A 



We can make the right-hand side (RHS) of ( fTSl l arbitrarily 
close to 1 by choosing arbitrarily small positive values for (. 
Thus we obtain (|4]i. 



3.2. Constant Nonzero Dark Current 

In this subsection we shall prove the first inequality in (O. To 
this end, we lower-bound on the mutual information /(Q*", W) 
for the input distribution ^ as follows: 

I{Q\W) 

= H{Y) - H{Y\X) 

OO 

= -Y^{(l-p)Vx{y)+prx+c(y)) 
y=o 

■\ogiil-p)rx{y)+pVx+ciy)) 

OO 

+ il-p)J2'P^iy)^ogVx{y) 

y=o 

OO 

+p^7'A+c(y)iog7'A+c(y) 

y=0 

|P«(.).0.((1-P)^+.) 
- (1 - P) f:v,<.y) log ((1 - p) +J>^^) 



OO 

= -p'^'Px+ciy) 

- (1 - p) £ r^iy) log (^(1 - p) + f 
= P E Pa+c (2/) log - E f (1_£)^ 

y=0 ^^^^ y=0 \ 

+ pP^Hm log ((1 -P)+ P^^) 

>0 V 



P ''aCb) y-^ 



- E ((1 - +P^A+c(y))p^^ 

/ -(C+A) (C+A)" X 

= PE^A+C(y)l0g .A,^' 

2/=0 \ ^ V'. / 

-ii-p)p±r,iy)'^ 



(1-p)pE^a+c(2/) VE 



e-(C+A)(C±A)^ 

c-A A" 



y=o 



y=o 



^pY,'Px+dy)(-C + v'^og(i + ^)) -(i-p)p 



y=o 



_p2 I y-e-(A+2C) 



E' 



A + 2C- 



2\ a 



2/! 



X 

OO OO 

KE^^+«(y) +pE^^+c(j')yiog ( 1 

i/=0 ■y=0 



= 1 

2 2 — 

\P + P — p e ^ 



= (C+A) 



: p(C + A) log ( 1 + ^ ) - K - P - ( - 1 ) . (17) 



For small enough we choose C = y A log ^ and p = ^ ^ 
= . By using (fTOl i and ( fTTb and letting 5 tend to zero 



\/Alog t 

we establish the lower bound in ^ 



4. THE UPPER BOUNDS 

In this section we shall sketch the proofs of the upper bounds 
on the asymptotic capacities of the Poisson channel. We shall 
use the duality bound IH which states that, for any distribu- 
tion R{-) on the output, the channel capacity satisfies 



C <snvE[D{W{-\X)\\R{-))] , 



(18) 



where the supremum is taken over all allowed input distribu- 
tions. We shall describe the choices of R{-) that lead to our 
upper bounds, but we shall omit the details. 

4.1. Dark Current Proportional to £ 

In this subsection we shall sketch the proof for (|5]l. To this 
end, as in 13], we introduce the Poisson channel with contin- 
uous output whose input x is the same as the original Poisson 
channel, and whose output is ?; e . The conditional den- 
sity W{-\-) is 

W{y\x)=V^+.A[y\)- (19) 

We denote the capacity of ( fT9l ) under Constraints (|2|l and ((Sj 
by C{\£,A). It is shown in |2| that 

C{\£,A) = C{X,£,A). 



Thus, to prove (|5]l, it suffices to prove 

<i, 0,00 . 

£[0 t log J 



(20) 



To this end, we choose the distribution ) on Y to be of 
density 

0<^<1 

where /? > is arbitrary, v £ (0, 1] and p G (0, 1) will 
be specified later, and r(-, •) denotes the Incomplete Gamma 
Function given by 

Applying ( fTSb on C(0,£,y^) with the above choice of fj^{-) 
in the place of R{-) and with the choice = ^ yields that, for 
every p S (0, 1) and (3 > 

C(0,f,^) <flogi + log-^ + | 

p I — p p 

+ £max|o,(ilog^ + log^ + i^)}. 

Choosing p = in the above inequality and letting E tend 
to zero yield (|5]|. 



4.2. Constant Nonzero Dark Current 

In this subsection we shall sketch the proof of the upper bound 
in ([8]). We choose the distribution R{-) on the output Y to be 

R{)^[^'^^^ ye{0,l,...,7V-l} 
\6{l~p)py-^, y£{N,N + l,...}, 

where N E Z+ and p E (0,1) are constants to be spec- 
ified later, and (5 is a normalizing factor given by <5 = 



y 



00 „-\\y 

--N ^ y\ 



We next apply ( fTsT i to upper-bound C{\,8 ,A). 



Calculation (with repeated applications of the Chernoff 
bound) yields 

< ( N log N + ^— + - log(27riV) + log 
- V 2 ' l-p 



£ 



N-VN - \ 



+ exp {N + NlogX- N log N) 



£ 



\£ 



1 + log - + log A , , 

P J V N-VN-X 

^ ^ ^ gAr-l-A+(A'-l)log A-(jV-l)log(A'-l) 



1 



A 



N - X 



O,log 



N - X 



(21) 



For small enough £, we choose N = [log j-J and let p E 
(0, 1) have any fixed value that does not depend on £. Apply- 
ing these choices to ( l2Tl i and taking the limit £ i yield the 
second inequality in (O. 
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